function diff = q_g_fun(m,param)

    diff = q_0_fun(min(m,param.m_u),param).*q_1_fun(min(m,param.m_u),param);

end